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Morphology mediates the interplay between the structure and electronic transport in atomically 
thin nanoribbons such as graphene as the relaxation of edge stresses occurs preferentially via out- 
of-plane deflections. In the case of end-supported suspended nanoribbons that we study here, past 
experiments and computations have identified a range of equilibrium morphologies, in particular 
for graphene flakes, yet a unified understanding of their relative stability remains elusive. Here, we 
employ atomic-scale simulations and a composite framework based on isotropic elastic plate theory 
to chart out the morphological stability space of suspended nanoribbons with respect to intrinsic 
(ribbon elasticity) and engineered (ribbon geometry) parameters, and the combination of edge and 
body actuation. The computations highlight a rich morphological shape space that can be naturally 
classified into two competing shapes, bending-like and twist-like, depending on the distribution of 
ripples across the interacting edges. The linearized elastic framework yields exact solutions for 
these rippled shapes. For compressive edge stresses, the body strain emerges as a key variable that 
controls their relative stability and in extreme cases stabilizes co-existing transverse ripples. Tensile 
edge stresses lead to dimples within the ribbon core that decay into the edges, a feature of obvious 
significance for stretchable nanoelectronics. The interplay between geometry and mechanics that 
we report should serve as a key input for quantifying the transport along these ribbons. 



INTRODUCTION 



The performance of nanoelectronic devices based on 
atomically thin films such as graphene depends critically 
on the interplay between geometry, structure and me- 
chanics. This is especially true in the case of nanorib- 
bons where the edge structure can fundamentally alter 
the overall response. As a classic example, the band 
gap in graphene nanoribbons (GNRs) is sensitive to edge 
type and ribbon width and in extreme cases determines 
the nature of the electronic transport, metallic or semi- 
conducting [THH]. The interplay with ribbon morphol- 
ogy serves as a crucial ingredient in quantifying these 
structure-property relations as the edge elasticity varies 
significantly with its structure [3 Ej. The edge shape 
also determines the nature and extent of edge function- 
alization [51 [TO]. Past studies on these ultra-thin rib- 
bons and sheets have revealed several interesting mor- 
phologies. In instances where the edge stress is compres- 
sive, the edges warp out of plane; Shenoy and cowork- 
ers used a combination of atomic-scale computations and 
scaling arguments to calculate the wavelength, ampli- 
tude and penetration width of such undulations in semi- 
infinite edges [UJ. Independent computations by Bets 
and Yakobson show that the rippling wavelength scales 
with the ratio of the edge stress to the flexural rigidity 
r e /D [12J . Below a critical width, the ripples transition 
into a spontaneous twist. The situation is expectedly 
different for tensile edge stresses, recently observed in 
reconstructed edges in graphene and intrinsic to bilayer 
(and possibly multi-layer) ribbons that reconstruct into 
partial edge tubules [T3J [14] . Here, the out-of-plane dis- 
placement occurs preferentially away from the edge such 
that the ribbon midsection curls as it ripples [T5] , 



Our focus here is on end-supported suspended nanorib- 
bons shown schematically in Fig.fTJi. Unlike freestanding 
nanoribbons, this architectural motif is a natural build- 
ing block for next-generation nanoelectronic devices and 
NEMS devices as it allows controlled yet scalable de- 
vice integration while minimizing deleterious substrate 
effects [IJ El [T51417J . The ribbons can be trimmed 
to shape before or after clamping them onto the end- 
supports (electrodes) [18-25], and the two scenarios re- 
sult in differing mechanical constraints on the ribbons, 
as detailed later in this article. Little is known regarding 
their morphological stability for several reasons. One, 
the edge stress-induced morphologies reported in past 
studies have been analyzed primarily to understand post- 
buckled shapes, via atomic-scale computations or scaling 
analyses with attendant simplifications [TTJ [J3]. How- 
ever, a detailed understanding of the pre-buckled shapes 
is necessary. These become tractable only in the (lin- 
ear) small amplitude limit and a naive approach would 
discount their utility as the post-buckling in these atom- 
ically thin sheets is expected to take place primarily 
through bending - the stretching is prohibitively expen- 
sive [26H29] . This is in stark contrast to recent compu- 
tations on edge morphologies of semi-infinite graphene 
sheets that show that stretching plays crucial, if not de- 
cisive role [11 . Furthermore, as we show below, the com- 
puted buckled shapes are quite sensitive to the initial 
perturbations, even more so for confined systems such 
as end-supported nanoribbons. Then, the stability of 
the precursory small amplitude deformations, where both 
stretching and bending can influence the stability, is the 
key to understanding the post-buckled shapes. Two, the 
end-supports as well as the edge-edge interactions limit 
the possible morphologies as they force to ribbon to be 
on average flat. More specifically, global buckling modes 
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such as (developable) twist or saddle shapes become un- 
tenable and edge stress accommodation takes place pri- 
marily through periodic ripples that also interact across 
the edges. Lastly, reconciliation of the shapes observed 
in experiments (for e.g. [T]d and [TJ;) or atomic-scale sim- 
ulations with elastic theories must factor in thermal ef- 
fects that naturally manifest as shape fluctuations. Sta- 
tistical theories of polymerized elastic membranes can be 
invoked to account for these effects via renormalized rib- 
bon elasticity [30] OS] ■ While we do not address this spe- 
cific issue here, it must be noted that the highly strain- 
sensitive bonding that stabilizes these flakes also results 
in strongly size-dependent and non-classical effects which 
can become significant at these scales [32] . 

In this article, we use stability analyses and compu- 
tations to explore the morphological stability space of 
nanoribbons as a function of both intrinsic and engi- 
neered parameters, i.e. ribbon geometry (width w and 
length I), material parameters (sheet and edge stiffness 
S and S e , and r e ) and deformation along the ribbon (x- 
)axis (uniaxial strain e xx ). Atomic-scale simulations of 
ribbons serve as inputs for identifying the possible equi- 
librium shapes. Although the computations focus exclu- 
sively on GNRs, we demonstrate the generality of these 
results via rigorous stability analyses based on isotropic 
elastic plate theory, which in turn allow us to develop sta- 
bility diagrams for a combination of geometric, material 
and processing parameters. 

PAST EXPERIMENTS AND COMPUTATIONS 

[l] showcases the periodic ripples that have been ob- 
served in suspended nanoribbons. The experimen- 
tal images ([T]d and [TJ;) show two differing morpholo- 
gies, bending-like rippling in multi-layered graphene rib- 
bons [33J and edge rippling in hollow BN ribbons [31]. 
The results of our atomic-scale simulations on graphene 
nanoribbons, summarized in [l];, reveal a considerably 
richer morphological space. To facilitate direct com- 
parison with past computational studies, we have per- 
formed these computations using a reactive bond order 
(AIREBO) potential [35] as implemented in the software 
package LAMMPS [35] (see Methods). For each combi- 
nation of parameters, the characteristics associated with 
the rippling are extracted from the relaxed shapes in the 
computations. The ratio of the width to wavelength A 
associated with the rippled shape is indicated in the fig- 
ure, expressed as a dimensionless wavenumber kw, where 
k = 2n/X. While we do not discuss the effect of rib- 
bon length I in detail in the remainder of the article, its 
main effect is to act as a constraint on the permissible 
wavenumbers in instances where n/l > k. 

Some of the shapes we observe are similar to those 
reported in past studies, both at zero and finite temper- 
ature [TTHT31 [37] . A key feature is that we observe two 



dominant morphological classes for each given set of pa- 
rameters: in-phase or symmetrical ripples (S), and out- 
of-phase or asymmetrical (AS) ripples, where the phase 
refers to the relative displacements of the edges. The spe- 
cific form of ripples is controlled by the sign of the edge 
stress and the aspect ratio. For compressive stresses and 
sufficiently large widths, the edges do not interact and 
the ribbon exhibits classical edge ripples (not shown) . As 
the width is reduced, the edge-edge interactions become 
important such that the ripples penetrate through the 
width and ribbon buckles either in- or out-of-phase [37] . 
At much smaller widths, the entire ribbon buckles with 
bending-like or twist-like undulations. Applied longitu- 
dinal strains e xx modify the midline-line morphology, as 
shown in the figure for symmetrically rippled ribbons. 
Tensile edge stresses force the midline to curl out-of- 
plane. However, since end-conditions require the midline 
to be flat, we observe dimples at and around the mid- 
line which decay into the edges. The dimples can also 
split asymmetrically at large widths, although the mode 
is relatively rare. At small widths, though, the ribbon 
morphology is either flat or symmetrically rippled. 

We note that some of these shapes, in particular sym- 
metric ripples in ribbons with compressed edges (e.g. 
filament-like buckling at small widths) are also observed 
in naturally occurring ribbons such as straight-edged long 
leaves j27j [38] [39] . While there are parallels with the 
rippled shapes analyzed here, the shapes in the materi- 
ally homogeneous natural systems are driven by inelastic 
and distributed growth strains unlike the highly local- 
ized and elastically stressed edges that characterize the 
nanoribbons considered here. Coupled with the differing 
end-conditions that arise in suspended nanoribbons, the 
stability of the expected shapes is qualitatively different. 

STABILITY ANALYSIS OF PERIODIC RIPPLES 

The diverse morphologies observed in the computa- 
tions are systematically analyzed using classical elastic 
plate theory. Consider a nanoribbon of thickness h and 
width w (h <C w) clamped to supports spaced apart 
by a length I. We assume that the ribbon is a lin- 
ear, isotropic elastic thin plate with elastic modulus E, 
Possion's ratio and bending and stretching stiffnesses 
D = Eh 3 /I2{l-v 2 ) and S = Eh, respectively. The edge 
stress r e in these atomically thin nanoribbons arises due 
to structural changes or reconstructions localized at the 
edges. In continuum limit, then, the edge can be approx- 
imated as an elastically stressed bounding spring with 
negligible bending stiffness [10] . This composite approxi- 
mation is similar in principle to the core-shell framework 
often invoked to describe elastic behavior or nanowires 
and thin films [1TH15] . For a system so structured strain 
compatibility and force equilibrium require that the elas- 
tic Hamiltonian that maps the initially flat ribbon to its 
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FIG. 1: (a) Schematic illustration of a suspended nanoribbon, simply supported or clamped depending on the end-supports. Ex- 
perimental observations of (b) rippling in multilayer graphene ribbons [33] and (c) edge rippling in hollow BN nanoribbons [34] . 
The images are reproduced with the authors' permissions, (d, e) Similar shapes observed in atomic-scale computations on 
I ~ 20 — 25 nm long, armchair-terminated graphene nanoribbons (a-GNRs) with (d) compressive and (e) tensile edge stresses. 
The critical wavenumber k*w is the ratio of the ribbon width to ripple wavelength. In all cases, both bending-like (symmet- 
ric, S) and twist-like (antisymmetric, AS) shapes are observed. The color indicates the scaled magnitude of the out-of-plane 
displacements as an indicator of the extent of edge-edge interactions. The additional effect of applied compressive and tensile 
strains for r e < is shown in bottom two plots in (d). The insets show details of the bent morphology at and around the 
ribbon mid-line. 



deformed state, R = (x, y, 0) — > R' = (x + u x , y + u y , £), 
satisfies the generalized Foppl-von Karman (F-vK) equa- 
tions [40] . The stability of the periodic ripples can be 
analyzed by assuming a sinusoidal variation in the out- 
of-plane deflection, 



C(x,y) = /(y)sinfcr. 



(1) 



The shape satisfies the boundary conditions along the 
simply supported sides of the graphene nanoribbon since 



C = and C,xx + vC 
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for x = and x = I. In 



the limit of negligible transverse and shear stresses, the 
classical F-vK equations simplify to a boundary value 
problem for the out-of-plane deflection (see Methods), 
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where e xx = e xx + e% x is a general uniaxial strain, ex- 
pressed as the sum of an intrinsic strain e xx , if present, 
and applied uniaxial strain e xx . Together with the 
boundary conditions, 



{[/, 



(f,yy~vk 2 f)\ ±f =0 

(2-^ 2 / ;5 ]T^ 2 /}|±f =0, 



(3) 



we completely specify the form of the deflection. 



The effect of the residual edge stress r e = t° + S e e xx 
that enters into the boundary value problem depends 
on the synthesis procedure. Here, is the unrelaxed 
edge stress. In cases where the ribbons are very long, 
or trimmed to shape following synthesis and then trans- 
ferred onto the end supports, the edge stresses result in a 
residual longitudinal force that must be borne by the rib- 
bon core. A simple force balance for the initially flat rib- 
bon yields the intrinsic body strain, 2r e +Swe xx — T ss 0, 
where T is the net longitudinal force. The body strain 
acts much like an imposed uniaxial strain over the entire 
ribbon core and therefore is analyzed in the context of 
cxtrinsically strained ribbons. An entirely different sce- 
nario occurs when the precursor flake is placed on the 
end-supports and then trimmed to shape. In this case, 
the edge-supports modify the net force balance such that 
the ribbon core is unstrained, e xx = 0. 
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FIG. 2: (a) Morphological stability diagram showing the scaled edge stress r e w/D as a function of the wave number kw for 
relaxed (e xx = 0, light shading) and intrinsically strained (e^ = —2t e /(Sw), dark shading)) nanoribbons with compressive 
edge stresses. Here as well as in the following figures the critical curves for bending-like and twist-like buckling are indicated 
by black and red lines, respectively, (inset) Magnified plot for small wavenumbers kw that shows the bifurcation between the 
two morphological classes for strained ribbons, (b) The ribbon profile obtained from the analytical solution (see SI) for relaxed 
ribbons plotted as the scaled deflection f(y)/f m ax versus scaled width yjw for the three different wavenumbers. For both 
classes, the rippling localizes to the edges with increasing wavenumber. Schematic illustrations of the shapes for some of the 
profiles are also shown. All plots are based on Poisson's ratio v — 0.17 corresponding to that for graphene [44| . 



COMPRESSIVE EDGE STRESSES 
Relaxed Ribbons 

We first analyze relaxed ribbons characterized by e xx — 
and r e < 0. The boundary conditions ( [3]) yield the 
critical point for the onset of buckling, conveniently ex- 
pressed as stability diagrams that relate the (dimension- 
less) edge stress r e w/D to the scaled wavenumber kw. 
[2£i shows these stability diagrams for both bending-like 
(S) and twist-like (AS) ripples in a ribbon with Pois- 
son's ratio corresponding to that for graphene, v — 0.17. 
The critical stress varies linearly both morphological 
classes at large widths (or equivalently short ribbons), 
t*w/Dk -[(l-v)(3 + v)](k*w)/2 for k*w > 1. The so- 
lution also yields the ribbon shape (SI, Eqs. S4 and S6), 
plotted in [2Jd as a shape function f(y)/f m ax fo r several 
representative wavenumbers. For small aspect ratios (e.g. 
k*w — 50), the ripples are uncorrelated and localized to 
the edges, as expected. At the critical point k*w s» 5, we 
see a bifurcation due to edge ripples that now begin to 
interact across the width via saddle-shapes morphologies 
with net negative Gaussian curvature, apparent in the 



schematic illustrations ui[2]d and also in the post-buckled 
shapes observed in computations ([TJi, with kw = 2.5). 

The asymptotic behavior for the limit kw <C 1 sheds 
light on the markedly different behavior for the two mor- 
phological classes. The bending-like ripples exhibit a 
quadratic dependence, t*w/D s» —(1 — v){l + u){k*w) 2 
while the critical edge stress for twist-like ripples is in- 
dependent of the wavenumber, t*w/D rs —4(1 — v). 
As an example, the almost flat (scaled) profile of a 
high aspect ratio ribbon with k*w = 0.5 is plotted 
in [2]d. The stability curve for the shape also follows 
from simple scalings based on elastic energies (calcu- 
lated per ripple wavelength A) associated with ribbon 
bending £(,, and ribbon stretching at the core and edge, 
£ s and E e s . The bending energy follows from the cur- 
vature tensor, Sj, ~ Dkw J Q (n xx + Ky y ) dx and £& ~ 

Dkw J Q A n xy dx for the bending-like and twist-like rip- 
ples, respectively. The stretching energies are related 
to longitudinal and edge strains, £ s ~ —Tkw J* e xx dx 

and £ s ~ —2r e kw J e xx dx. Ignoring the weak varia- 
tion in out-of-plane deflection across the width, / ~ S x 
and / ~ 2y/w5 x for the bending-like and twist-like 
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FIG. 3: (a-b) Comparison of the analytically predicted stability of of bending-like and twist-like shapes and the buckled 
shapes observed in atomic-scale simulations for (c) relaxed and (d) intrinsically strained a-GNRs, now plotted as k(w). The 
representative atomic-configurations for symmetrically rippled ribbons with widths w = 0.72 nm, w — 1.93 nm and w — 5.12 nm 
are also shown (inset). 



ripples, respectively. Then, k xx ~ —k 2 S x smkx and 
K yy ~ and the dominant contributions for bending- 
like ripples scale as £b ~ Dk 4 S 2 w and £| ~ —T e k 2 S 2 . 
Taken together, they yield the quadratic dependence, 
—T*w/D~(k*w) 2 . In the case of twist-like ripples, the 
stretching energy remains unchanged. The bending en- 
ergy is due to K xy ~ k5 x /w cos kx and is relatively larger, 
£b ~ Dk 2 /w. Equating the two yields the critical edge 
stress, —t*w/D~1. 

The stability diagram yields insight into the shapes ob- 
served in atomic-scale simulations on GNRs. For unre- 
constructed armchair and zigzag terminations, the com- 
pressive edge stresses range from t® = — 10.5eV/nm 
to — 20.5eV/nm while the edge stiffness S e varies from 
113 to 147eV/nm H]. The low bending stiffness of 
these atomically thin sheets (D ~ 1.5 eV) yields a scaled 
edge stress r e w/D ~ 10 for nanometer- wide GNRs. The 
critical wavenumber is in the vicinity of the bifurcation 
point where both bending- and twist-like undulations are 
possible, although the former are energetically favored. 
This is corroborated by the computed morphologies for 
a-GNRs with varying widths, shown in[ljl-e and [3^ (in- 
set). Both morphological classes are observed depending 
on the form of the perturbation. [3^, shows the reduced 
stability diagram, k vs. w predicted by our analysis 
and that extracted from the computations. The over- 
all trend is well-described by our stability diagram al- 
though there arc quantitative deviations in the simulated 
shapes. They are likely due to the continuum approxi- 
mation of an atomic-scale system and the fact that we 
ignore the relaxation of the edge stress due to out-of- 
plane displacements [8]. Furthermore, a more detailed 
analysis of the atomic-configurations reveals that at in- 
termediate to large widths, the two morphological classes 
can also co-exist, albeit with differing ripple wavelengths. 
This can be clearly seen in the atomic-configuration for 



the w = 1.93 nm a-GNR shown in Fig. [3^, (inset). We see 
an out-phase component in the symmetric ripples away 
from the ribbon ends, indicative of a twist-like buckled 
mode with a much longer wavelength that allows the rib- 
bon to further relax the residual edge stress. At very 
small widths, the co-existence shapes lead to locally flat 
morphology quite like the one shown for w — 0.72 nm. 

Strained Ribbons 

The eigenvalue solution can be written as 

$(kw, T e w/ D, e xx Sw 2 /D) = 0, (4) 

where e xx Sw 2 /D is the additional (scaled) strain. As 
a starting point, we explore the effect of intrinsic body 
strains that arise naturally in long nanoribbons, e xx = 
e xx . The compressive edge stresses require that e xx — 
~2r e /S > 0. The modified stability stability diagram, 
also shown in[2^,, is based on a lengthy analytical solution 
which is not shown for brevity (see SI) . It is immediately 
clear that the flat phase is more stable as the edge-edge 
interactions can be absorbed for relatively wider ribbons 
due to the mediating strained core. In effect, the relax- 
ation of the edge ripples shifts to smaller wavenumbers 
and larger edge stresses as it now occurs with respect to 
an already stretched ribbon core. The bifurcation into 
the two morphological classes shifts accordingly (inset). 
Of importance is the dramatic effect of the intrinsic strain 
on the relative stability of the bending- and twist-like rip- 
ples for small aspect ratio ribbons, kw <C 1. Although 
the critical point for bending-like ripples is still lower 
(inset), we no longer see the quadratic dependence in the 
critical point for bending-like ripples. Rather, quite like 
the twist-like rippling in relaxed ribbons, both morpho- 
logical classes result in rippled states that are determined 
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entirely by the (larger) critical stresses, t*w/D, that are 
in turn independent of the wavenumber. The fact that 
twist-like ripples compete effectively is not surprising as 
they preserve the midline length and therefore are clearly 
efficient in accommodating the tensile strain within the 
core. In the specific case of a-GNRs, the predicted stabil- 
ity space k(w) as well as the results of the computations 
are plotted in[3]D. As in the unstrained case, the analyti- 
cal framework captures the trends although the quantita- 
tive agreement becomes noticeably poor at small widths 
due to effects mentioned earlier. 

The interplay between bending-like ripples and intrin- 
sic body strains can be understood from the underlying 
energetics at these intermediate widths. Approximating 
the saddle morphology as f(y) s» S x + [1 — cos(ny /w)]S y 
with a transverse amplitude constant much smaller than 
that along the ribbon, \S y \ <C 5 X , the bending energy 
due to the two curvatures k xx ~ —k 2 5 2 sinfcx 
S y /w 2 smkx scales as £i> ~ Dw(k A 5 2 + dy/w 4 ). The 
stretching energy has contributions from both the core 
and the edge, £ s ~ i~ e k 2 d 2 and £f ~ —T e k 2 (5 x + 8 y ) 2 . 
Then, the leading order term for the total energy for 
stretching is £ s + £% ~ —T e k 2 S x 6 y . The relation between 
the two amplitude constants follows from the zero net 
moment across the ends, which reduces to 

— Dwk xx + Td x sin kx + 2r e (5 x + 5 y ) sin kx ~ 0, (5) 

where the longitudinal force T = — 2r e . Simplifying, 
we get the ratio S y /S x ~ — (kw) 2 /(T e w/D) which is 
also related to the intrinsic (scaled) strain, S y /S x ~ 
(kw) 2 /(e xx Sw 2 /D). Substituting in the balance, £b = 
£ s + £g yields the observed dependence, t*w/D ~ — 1. 

The stability space for extrinsically strained nanorib- 
bons is shown in [4^, and [4Jd for bending-like and twist- 
like rippling, respectively. Note that the compressive 
strains e xx < result in a different class of solutions 
as detailed in SI. The contours in the plots correspond to 
critical (scaled) edge stresses as a function of (scaled) ap- 
plied strain and wavenumber. The iso-(edge)stress con- 
tour for classical Euler buckling (r* = 0) and iso-strain 
plots for intrinsically strained nanoribbons (e xx — e xx ) 
are also plotted. For small strains and large kw ^> I 
characterized by non-interacting rippled edges (5 y ~ 0), 
the critical strain decreases linearly with wavenumber. 
Here, the bending energy due to the rippled edges is con- 
strained to a penetration width wi ~ 1/k and there- 
fore scales as £b ~ Dk 4 wi S 2 ~ Dk 3 S 2 . Both the ribbon 
core and the edge contribute to the stretching energy, 
which scales as £ s + £% ~ — {e xx Sw + 2T e )k 2 5 2 . Equating 
the two yields the linear dependence, i.e. e xx Sw 2 /D ~ 
— (k*w + T*w/D). At small wavenumbers kw < 1, on 
the other hand, the critical applied strain is proportional 
to the edge stress and is independent of the wavenum- 
ber for both morphological classes. This follows from 
the analytical solutions as well as scaling analyses for in- 
trinsically strained ribbons where the interplay between 



transverse and longitudinal curvatures results in a critical 
edge stress t*w/D that is independent of the wavenum- 
ber ([2]); the critical strain varies similarly as evident from 
the stability diagrams. [4"5] 

The transition behavior at intermediate wavenumbers 
kw ~ 1 is characterized by peaks in the critical strains 
for ribbons subject to tensile and compressive strains. 
The peaks are apparent in the plot for bending-like rip- 
ples and appear at strains larger than the range shown 
in the plot for twist-like ripples. |4j: shows the iso-stress 
stability curves for the specific case of w = 10 nm wide a- 
GNRs (edge stress r e = — 10.5eV/nm). The wavenumber 
is plotted on a log-scale to highlight the relative stabil- 
ity of the two morphological classes at small wavenum- 
bers. At applied strains less than the peak strain we see 
a bifurcation into an additional morphology with smaller 
wavenumbers. The double-well profiles associated with 
this sub-class (k*w = 7.5) are shown in[4jl for a-GNRs 
subject to a uniaxial strain e xx — 1%. The co-stable mor- 
phology with larger critical wavenumber exhibits ripples 
localized to the edges (k*w = 37), as expected. At much 
lower strains, we recover the saddle-shapes analyzed ear- 
lier. In the case of a-GNRs, the transition occurs at 
strains below w 0.3%. As confirmation, the bending- 
like rippled morphology for e xx — 0.15% and k*w = 0.5 
is shown in [4ji. 

The additional double- well shapes are also observed in 
our computations. In fact, transversely rippled center- 
line morphology is clearly visible in Fig. [TJi for bending- 
like ripples in ribbons subject to tensile strains, e xx = 1%. 
It is interesting to note that these profiles are more sus- 
ceptible to twist-like rippling as the ribbon midline is 
largely unmodified compared to bending-like ripples. In 
fact, for a range of applied strains (0.1%-0.3% for the a- 
GNRs), the double- well shaped rippling can only occur 
in a twist-like fashion (Kb). Similar peaks also occur in 
critical strains for ribbons with tensile edge stresses and 
subject to compressive applied strains. The behavior is 
consistent with classical Euler buckling of ribbons appro- 
priately modified due to the edge stress; the more general 
morphology of ribbons subject to tensile edge stresses is 
discussed in the next section. 



STRAINED RIBBONS WITH TENSILE EDGE 
STRESS 

Tensile edge stresses are usually driven by reconstruc- 
tions and therefore can vary significantly. Recent calcu- 
lations indicate that this is indeed the case in graphenc 
nanoribbons. a-GNRs with 5 — 6 reconstructions re- 
sult in large stresses, r e = 24eV/nm while those in 
the 5 — 7 reconstructed z-GNRs are almost negligible, 
t p = 0.02eV/nm. The periodic ripples in the sus- 
pended nanoribbons can again be categorized broadly 
into bending-like and twist-like, and both classes have 
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FIG. 4: (a-b) Iso-(edge)stress critical contours for (a) bending-like and (b) twist-like rippling of strained ribbons with compres- 
sive and tensile edge stresses. The contour plot associated with classical Euler buckling (r e = 0) is also plotted. The solid white 
lines are iso-strain plots for intrinsically strained nanoribbons summarized in [5] and [Sp. (c) The critical curves for strained 
unreconstructed a-GNRs with compressive edge stress r e = 10.5eV/nm. (d) Ribbon profiles f(y)/f m ax for some of the critical 
points below (e xx = 0.15%, k*w — 0.5) as well as in the vicinity of the peak strain (e xx = 1%, k*w — 7.5, 37). 



been observed in simulations ^f)- The solutions for the 
rippled shapes in the small amplitude limit follow directly 
from T]|3 and are again detailed in SI (Eq. S7) 

We restrict our analysis to compressively strained rib- 
bons as the morphology of unstrained or stretched rib- 
bons is trivial - they are always flat since the edge stress 
is completely absorbed by the relevant end conditions. 
The reasoning is consistent with our analysis in that the 
critical contours for rippling shown in Fig. [4^,-b always lie 
below the e xx = line; the rippling requires a compres- 
sive stress. For intrinsically strained nanoribbons, the 
iso-strain contour line for bending-like ripples is plotted 
in Fig. (solid white line below r e > 0) while the de- 
tailed stability diagram r e w/D vs. kw and the associated 
ribbon profiles are plotted in [5j At large wavenumbers 
the edges do not interact and past the critical point the 
core develops dimples around the midline (positive Gaus- 
sian curvature) to accommodate the compressive strain. 
The limit corresponds to Euler buckling of the ribbon 
core with a scaled critical edge stress that varies quadrat - 
ically with the critical wavenumber, as expected. The 
edge interaction at smaller wavenumbers enhances the 



relative extent of the dimples, i.e. the transverse curva- 
ture increases. In the limiting case kw -C 1, the critical 
wavenumber is independent of the edge stress. The trend 
is qualitatively similar to that in ribbons with compres- 
sive edge stress (|2^) for the simple reason that in both 
cases the bending energy is modified by transverse cur- 
vatures, albeit opposite in sign. 

The twist-like rippling in these ribbons consists of 
dimples that alternate across the ribbon midline, cor- 
responding to an axial-torsional buckling mode of the 
ribbon core. At small wavenumbers (width), the criti- 
cal edge stress is expectedly larger due to the significant 
transverse bending (Gaussian) energy associated with the 
morphology. Increasing the wavenumber reduces this en- 
ergy and the critical stress decreases quadratically. In the 
large wavenumber limit, the edges no longer interact and 
the critical curve is identical to that for bending-like rip- 
ples, i.e. the critical stress increases quadratically with 
wavenumber. At intermediate wavenumbers, therefore, 
we see a minimum in the critical edge stress. Note the 
ribbon profiles are relatively insensitive to the wavenum- 
ber. The scaled edge stress for nanometer wide a-GNRs 



fall within the proximity of the minimum. Although the 
bending-like ripples are favored, the energetic difference 
is small in that we see twist-like ripples in some of the 
computed morphologies (fife). Finally, the stability with 
respect to extrinsic compressive strains is shown in [4]a,- 
b. Larger compressive strains always favor bending like 
ripples since the rippled mid-line is more efficient in ab- 
sorbing the imposed strain. Of note is a peak critical 
strain for bending-like ripples that reflects the transition 
into a higher order Euler buckling mode (in this case, an 
S-shaped profile) at larger wavenumbers. In fact, there 
are several such peaks that emerge at larger compressive 
strains and wavenumbers as the higher buckling modes 
become viable. 
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ductance fluctuations to band gaps. Clearly, the inter- 
play between ribbon structure, geometry and externally 
applied strains results in a morphological stability space 
that is considerably richer than the classical edge rip- 
ples that are normally associated with these crystalline 
nanoribbons. For example, tensile edge stresses tend to 
destabilize the ribbon core while forcing the edges to be 
flat, and the dimpled morphology affects both equilib- 
rium and transport properties. We note that the excel- 
lent agreement in the trends observed in the predicted 
morphologies and those observed in atomic-scale simu- 
lations in graphene nanoribbons is limited to low tem- 
peratures as we have ignored the effect of thermal fluc- 
tuations. This is to be expected as the phonon fluctua- 
tions that couple the in-plane and out-of-plane displace- 
ments renormalize the ribbon stiffnesses such they be- 
come inherently size dependent. The stability diagrams 
extracted here are still applicable in that the dominant 
effect of these fluctuations is to change the scaled edge 
stresses and strains. In the case of nanoribbons, though, 
there is an inherent anisotropy due large differences in 
the thermal fluctuations along the transverse and lon- 
gitudinal directions. As an interesting example, in sys- 
tems where the base bending rigidity increases, the twist- 
like dimples can stabilize at high enough temperatures as 
the rigidity is enhanced preferentially along the ribbon, 
D x > D y . In that statistical understanding of the 

effect of the resultant shape fluctuations on the ribbon 
properties becomes necessary. 
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METHODS 



Atomistic Simulations 



FIG. 5: (a) Stability diagram and (b) ribbon profiles as in [2] 
but for intrinsically strained nanoribbons with tensile edge 
stresses. For clarity, the dimpled morphologies predicted by 
our analysis for k*w = 5 and k*w = 50 are shown schemati- 
cally. 



CONCLUDING REMARKS 

The final ribbon morphology has ramifications for elec- 
tronic properties that are strongly coupled to the ripples 
and the associated strain distribution, ranging from con- 



The computations are performed on graphene nanorib- 
bons with a fixed length, I = 21 nm. For extrinsically 
strained ribbons, the strain was applied to the fixed ends 
of the ribbons. The width was systematically varied in 
the range w = 1.2 — 10 nm. The effect of compressive 
and tensile edge stresses was studied using pristing and 
reconstructed graphene edges, for armchair and zig-zag 
edge structures. In each case, the nanoribbons were per- 
turbed by a combination of random and sinusoidal dis- 
placements and then relaxed to their equilibrium shape 
using a conjugate gradient algorithm with an energy tol- 
erance of 10 -10 eV. 
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Composite elastic framework 

The nanoribbon is approximated as an isotropic elas- 
tic thin plate with thickness h, width w, and length I 
(/i«i«Cf). Our composite framework consists of edge, 
modeled as a stretched or compressed elastic string that 
is glued to the ribbon. The edge stress T e is the main 
material parameter associated with the string, it has no 
bending stiffness. Additionally, the edge is assumed to 
be sharp such that that the width of the ribbon is in fact 
that of the ribbon core. 

Governing Equations 

The governing equations follow from standard Fv-K 
equations in the small deflection limit, 

<1",,„,.., V 4 $ = S{C,xxC,yy ~ (%y)- (6) 

Here $ is the Airy function associated with the in-plane 
stresses and the operator \7 4 A — A tXXXX + 2A :Xxyy + 
A^yyy. For long ribbons, end effects can be safely ig- 
nored; the ribbon has negligible transverse and shear 
stresses, and the governing equation simplifies to 

.XXX X + 2C ,xxyy + Cyyyy) — TC,, XX 

where T = Se xx is the net longitudinal force. Setting 
our origin on the midline, the force balance modified by 
the edge stresses and free torques specify the boundary 
conditions, 

{t,vv + Kxx}\±$ = 0, 

{c,„ + (2-WC,..±gc}l ± f =o. 

Substituting the sinusoidal deflection £(x, y) = 
f(y) sin kx in the equations yields the eigenvalue 
problem for the ribbon shape in terms of the ribbon 
strain e xx , Eqs. 2 and Eq. 3. 

Analytical Solutions 

The parameter sets considered here arc yield fully ana- 
lytical solutions for the ribbon morphology. In particular, 
the form of the solutions differs for relaxed and strained 
ribbons, as detailed in Supplementary Information (SI). 
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For each parameter set (r e , e xx ), our composite framework yields fully analytical expressions for the critical ribbon 
shapes. The solutions for the parameter sets analyzed in the text are detailed below. 



ZERO STRAIN e xx = 



The governing equation for relaxed suspended nanoribbons reduces to 



f,yyyy 2fc 2 / iTO + k 4 f — 0. (SI) 
The general solution for both bending-like and twist-like ripples is of the general form 

f=Yl A n(ky) n cosh ky + B n (ky) n smhky, (S2) 

71=0,1 

where the set of constants (A ,Bi) and (Ai,B ) correspond to bending-like (symmetric, S) and twist-like (anti- 
symmetric, AS) ripples, respectively. The constants can be evaluated using the boundary conditions (Eq. 5). The 
solution also yields the critical buckling point and the associated shape function. For bending-like ripples, we get 

t*w _ (l-v)k*w[{Z + v)smh.k*w-{l-v)k*w\ 

~D~ ~ 2(coshk*w + 1) ^ ' 

At small wavenumbers, kw < 1, t*w/D = -(1 - v)(l + v)(k*w) 2 and for kw > 1, t*w/D = -(1 - v)(3 + v)k*w/2. 
The associated shape function is, 

(1 - v)kwta,n\\ (kw/2) + 4 / ,n 
/ = A{ ± 2(1 -z/) C ° V ~ V V ^ ^ ^ 

where A is an arbitrary constant. The corresponding relations for twist-like ripples are 

t*w _ (1 - v)k*w[(3 + v) sinh k*w + (1 - v)k*w] 
~D~ ' 2(coshfc*w - 1) 



(S5) 



Again, for kw < 1, t*w/D = -4(l-v) and for kw > 1, t*w/D = -(l-v)(3 + v)k*w/2. The solution for the ribbon 
profile is 

a r (1 - ^)fcwcoth(fcw/2) + 4 . , , , , , , 
/ = A{ K - ' 2{1 _ v) sinh ky - ky cosh ky} (S6) 



TENSILE STRAIN e xx > 

The general solution for tensile strains e xx > is 

/ = A coshpy cos qy + B Q coshpy sin qy + Ai sinhpy sin qy + B\ sinh py cos qy. (S7) 

Here p = k^J{^l + kf/W + l)/2 and q = fcy / ( v /l + kJJW - l)/2, where k r is defined as k 2 r = T/D. As a check, 
the solution for e xx — 0, Eq. 6, can be obtained in the limit fc r <Cfc and by letting coshpy— > cosh ky, sinh py-^ sinh ky, 
cos qy— >1 and singz— >ky. 
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COMPRESSIVE STRAIN, e xx < 

In the case of compressive strains, the general solution for k r < k is 

/ = Aq coshpy + B\ sinh qy + Ai cosh qy + £?i sinh qy (S8) 
Here A:, 2 = —T/D, p = kyjl + k r /k, q = ky/l — k r /k. For small wavenumbers, the modified solution takes the form 

/ = A n coshpy + _B sinh qy + A\ cos qy + B\ sin qy (S9) 
where p = k^/k T jk + 1, g = k^/k r /k — 1. For special case k = k r , 

f = Aq cosh \/2fcy + _B sinh V2fcy + Ai + Siy. 

INTRINSICALLY STRAINED RIBBONS, e xx = -2r e /S 

The general eigenvalue solution for generally strained ribbons can be written as a function of scaled wave number 
kw, scaled edge stress r e w/D and scaled strain Tw 2 /D, 

^ = ^{k*w,T*w 2 /D) (S10) 

For freestanding nanoribbons, we can obtain the solution by substituting e xx = —2r e /Sw into Eq. S10. The results, 
t*w/D = $>(k*w), are shown graphically as white line in Figs. 3a and 3b of main text. Due to the complicated 
mathematical formulae for ^ {k* w ,T* w 2 / D) and $(k*w), t*w/D can be only evaluated numerically. 



